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The Goos-Hanchen (GH) shift and the Imbert-Fedorov (IF) shift are optical phenomena which 
describe the longitudinal and transverse lateral shifts at the reflection interface, respectively. Here, 
we report the GH and IF shifts in Weyl semimetals (WSMs)—a promising material harboring low 
energy Weyl fermions, a massless fermionic cousin of photons. Our results show that GH shift in 
WSMs is valley-independent which is analogous to that discovered in a 2D relativistic material— 
graphene. However, the IF shift has never been explored in non-optical systems, and here we 
show that it is valley-dependent. Furthermore, we find that the IF shift actually originates from 
the topological effect of the system. Experimentally, the topological IF shift can be utilized to 
characterize the Weyl semimetals, design valleytronic devices of high efficiency, and measure the 
Berry curvature. 


Introduction. —Modern quantum physics originates 
from understanding the wave-particle duality of all parti¬ 
cles, among which photon is the first one being discovered 
with such duality. Classically, the physics of a beam of 
light being reflected at an interface is governed by geo¬ 
metric optics law, where the photons are treated as clas¬ 
sical particles. In contrast, when considering the wave 
nature of photons, spatial shifts at the interface appear 
as longitudinal shift in the incident plan^IH^, or trans¬ 
verse shift normal to the incident plan^^K^, which are 
known as the Goos-Hanchen (GH) effect and the Imbert- 
Fedorov (IF) effect, respectively. Due to all particles pos¬ 
sessing the wave-particle duality, the spatial shifts are 
also expected for other particles. For example, GH effect 
has been shown to exist in the systems of electron^^^^^, 
neutron^^, atomJ^, etc. Particularly for the 2D mass¬ 
less Dirac fermions in graphene systems, the GH shift can 
be manipulated froni positive to negative by tuning an 
external electric fielcP^. However, the IF effect has not 
been studied in non-optical systems. 

Similar to photons, Weyl particles are also massless. 
But different from photons, Weyl particles are spin 1/2 
chiral fermions and described by the Weyl equation. Re¬ 
cently, Weyl semimetals (WSMs) have been proposed 
as promising systems embedding Weyl fermions, gen¬ 
erating intensive interestPH^. In WSMs, the Weyl 
nodes always exist in pairs with opposite chir alities , and 
each Weyl node corresponds to a valley inded^^^. Sev¬ 
eral candidates are suggested to be WSMs including py- 
rochlore irradiated, to^ogical insulator and normal in¬ 
sulator heterostructured, staggered flux states in cold 
atom systemd^, and photonic crystals based on double- 
gyroid structured^. Despite of various theoretical pre¬ 
diction of WSMs, the experimental realization of WSMs 
remain a challenge. This failure cannot be ascri bed to 
the impediment of material growth techniqud^^, but to 
the scarcity of the method of direct experimental iden¬ 
tification of WSMd. Due to the topological proper¬ 
ties of the Weyl fermions, the WSMs may possess ex¬ 
otic wave-packet dynamics, which indicates a new route 


to characteriziM WSMs and potential applications in 
valleytronicdM^. 

In this paper, we report the GH effect and IF effect 
for 3D Weyl fermions in WSMs. By using wave packet 
method, we derive analytic results for the spatial shift 
of the GH effect and IF effect. Our results show that 
the GH shift is valley-independent. By contrast, the 
IF shift is valley-dependent [see Fig. [^, which give rise 
to the valley-dependent anomalous velocities in the sys¬ 
tem. Due to the IF shift being perpendicular to the in¬ 
cident plane, it is a generic 3D effect and could never 
appear in a 2D material, e.g. the graphend^. Further¬ 
more, we demonstrate that the IF shift originates from 
the topological effect of the system, namely, the Berry 
curvature of the system. Remarkably, the consequence 
of the valley-dependent anomalous velocity is significant 
enough to be detected experimentally. Finally, we discuss 
three applications of the valley-dependent IF shift: (i) ef¬ 
fectively characterizing the WSMs; (ii) directly detecting 
the Berry curvature; (iii) efficiently inducing valley cur¬ 
rent with a high polarization rate. 



FIG. I: Illustration of GH effect and IF effect in WSMs. The 
orange arrow represents the incident wave packet which in¬ 
clude Weyl fermions from two valleys, whereas the green (red) 
arrow represents the reflected wave packet of Weyl fermions 
only from Valley A (B). Agh denotes the GH shift, and 
stands for the IF shift of valley A (B). This figure only shows 
the positive GH shift case. 

Quantum GH and IF effects in WSMs .—The Hamil- 
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tonian of the WSMs system is 


X; Vi Pi CTi {x ^ 0) 

V'-Piai ^V{x) {x> 0) ’ 


where k = {hVykyY + {hv'^kzY — {E — Vy/hv'^, /3 = 

7 = v'^lvz, a = {vyhy/v^k^jc), and 0 = 

tan“^ (ky/kx)‘ C = sgn[vxVyVz\ is the chirality of the 
valley in WSMs. 


where Vi {v[) is velocity parameter for region x ^ 0 
{x > 0); Pi are momentum operators, V{x) is poten¬ 
tial, and (7i stand for Pauli matrices. To guarantee 
the model Hamiltonian Hermitian, we let = v'^ in 
this paper. This Hamiltonian indicates an interface lo¬ 
cated at X = 0 in the WSMs[see Fig. I^a)]. We con¬ 
sider a beam of Weyl fermions incident from the re¬ 
gion X < 0 modeled by a Gaussian wave packet as 

V’f (r) = ST^dkydk^ f{ky-ky)f{k^-k;)^^'^{\^,Y), 

where /(fcg — kg) = g-(*s-fea) /2A^^ 

Gaussian distribution functions of width peaked at 
the mean wave vector {kx^ky^kz) with s = y^z. Note 
that none of our results depend on the shape of the wave 
packet. Here, r) is the incident wave function, 

which is a solution of Weyl equation, i.e., = Epj'^^ 

for region x < 0: 


r) 


1 f 

yi + r?2 V T] 


ikxX+ikyy+ikzZ 

^ 5 


( 2 ) 


where kx = ^/E‘^ — {^VykyY — {kVzkzY/hvx^ cx = 

ta,n-^{vyky/v:ck:g), and rj = ^ E+hllkl ■ Analogously, 
the reflected wave packet can be written as ' 0 ^^(r) = 

!Tooj-oo^^ydkzf{ky - ky)f{kz - k;,)'ip^^(k,r), where 
is the reflected wave function, p;^^ can be ob¬ 
tained from the incident wave Eq. § by the substitution 
kx ^ —kx^ a TT — a and multiplication with the reflec¬ 
tion amplitude r = \r\e^^^. The integrals of pj^^p and 
give the center of the wave packets, and therefore we can 
obtain the spatial shifts in y, z directionJ^ 

f ~ ~'^^r(ky, kz) =F '^^^{ky, kz) 

\^^± = -i^MkyVkz)T^a{ky,h) • ^ ^ 


The spatial shifts for Weyl fermions are defined as the 
average shifts of the two spinor components: = 

(A^G)^^ 2 a^G))/(i^^ 2 )^ InEq.(3), aiky^kz) represents 
the incident angle, and p)r is the phase of the reflection 
coefficient, which can be obtained by matching the wave 
function at x = 0. If the incident wave packet is confined 
in x-y plane, the in-plane shift A^ and out-of-plane shift 
A^ of wave packet corresponds to the GH shift /S^qh and 
IE shift /^iF' 


^GH 

Aif 


Vy (1 + /3 sin^ - ^) 

P sin a cos a n 

—C I I ^ ^ ~ V ~ 7 ) 

Vy E idXiO 1 -^( 1 -^) 


( 4 ) 

( 5 ) 



FIG. 2: (a) Schematic of p-n junction with an interface at a: = 
0. Vi {v'i) is the velocity in the left (right) side of the interface. 
E is the Fermi energy, and V is the potential difference of the 
junction, (b) GH shift Agh and IF shift Ap^^ of valley A (B) 
versus incident angle 0. Total reflection happens at ^ ^ 30° 
(the yellow color shadow region). \p = hvi/E stands for the 
wave length of Weyl fermions, (c) Potential dependence of 
the GH shift for incident angle 0 = 30° (solid), 45° (dashed), 
and 60° (dotted), (d) Potential dependence of the IF shift for 
incident angle 6 — 30° (solid), 45° (dashed), and 60° (dotted). 


We take the incident Fermi energy E = 100 meV, po¬ 
tential V = 150 meV, and velocities Vi = v[ = 10^ m/s, 
where i = x ,y, z. In this case, /3 = 7 = 1. Note that 
the two valleys (A/B) of a WSMs have opposite chiral¬ 
ities. Fig. I^b) shows the spatial shifts versus the inci¬ 
dent angle 6 >, where both the GH and IF shifts are odd 
functions of incident angle, which is consistent with sym¬ 
metry analysis. Fig. [^c) shows the valley independence 
and potential dependence of GH shift Aqh^ which can 
be tuned from positive to negative by external field V. 
This feature is analogous to the GH shift in graphen^^. 
Fig. [^d) illustrates that the IF shift is independent of 
potential, but depends on valley index. The IF shift can 
be utilized to manipulate the valley degree of freedom. 

Topological origin of the IF effect. —Based on the semi- 
classical dynamics of wave packet, we show that the IF 
shift is closely related to the Berry curvature of the sys¬ 
tem. Let us assume the velocity Vi = v[ (i = x, ^, z) in 
WSMs system. In this case, = 7 = 1, and Eq.([^ 
reduces to Ajf = tan e • Gommonly, the 

Weyl node can be regarded as a magnetic monopole in 
/c-spac^^, and thereby generates an effective magnetic 
field in /c-space. The Berry curvature of Hamiltonian 
Eq.([^ in region {x < 0 ) is for con¬ 

duction band and valence bancPl^ respectively. The semi- 
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FIG. 3: (a) Schematic of the valley splitter for Weyl fermions. 
Region I, II, and III are three WSM layers with different ve¬ 
locities v^, and (b) Illustration of wave packet 

trajectory (black arrow) of Weyl fermions in region II. The 
electrons are injected by TEM tip with incident polar angle 
0 and azimuthal angle cj). The parameter 0 is fixed 50° in our 
calculation. In ordinary case, the trajectory of wave packet 
stays in ABMN plane. However, due to the valley-dependent 
IF shift, the trajectory of Weyl fermions from valley A (B) 
would shift to plane ABRS (ABPQ). (c) 0 dependence of the 
total anomalous velocities, (d) The propagating plane shift 
for two valleys. After considering IF shift, SN (QN) is the 
final position shift for valley A (B). 


classical equation of motion (EOM)P^ of wave packet 
~ ^ ^ where r and k are the center 

positions of the wave packet in phase space. We assume 
that the incident wave packet locates in the conduction 
band, and consider the incident wave packet in x-y plane 
with kz = 0. The variation of kx by potential V and 
nonzero Berry curvature Qy leads to the IF shift Ajp in 
2 ;-direction: 


Aif = 


r 

Jk^ 


dkx — ( ) 


E tan 0 ’ 


(6) 


where tan^ = ky/kx- Remarkably, Eq. ® and Eq.@ 
completely coincide in the case of /3 = 7 . This co¬ 
incidence does not depend on linear dispersion of the 
systenP^. The consistency between Berry curvature cal¬ 
culations and wave packet results strongly support that 
the IF shift is mainly a topological effect. The IF shift 
in WSMs is quite different from that in optical systems. 
Particularly, Aif reaches the maximum in WSMs [see 
Fig. 2 (b)J comparing the zero value in optical systems at 
Q — 0 ®^. This is because the conservation of ky^ kz guar¬ 
antee the Weyl Fermions stay in the same valley during 
the reflection processes in WSMs. In contrast, the polar¬ 
ization of the photons change during the reflection pro¬ 
cesses, which will severely influence Ajf- However,the 


semiclassical equation cannot be generally applied to all 
cases of reflection process because of the breakdown of 
adiabatic approximation for some systems. 

Anomalous velocities induced by IF effect. —We 
consider a well collimated beam of Weyl fermions 
propagating in the middle layer (Region II) of a sand¬ 
wich structure, which is constructed by three layers of 
WSMs[see Fig. |^a)]. The applied electric potential pro¬ 
file is shown below the sandwich structure. The height 
of the structure is h and the width of the region II is d. 
There are both GH and IF shifts at the two interfaces 
(x = ±|). In order to observe the valley-dependent 
IF shift, the mirror symmetry about x = 0 plane need 
to be brokerP^. Thus, we consider that the z direction 
velocities in the three regions are different, with vl = 

= V, and = vr; whereas the x and y directions 
velocities are still identical 

Without considering the spatial shifts at the interfaces, 
the normal velocities in region II in y and z directions 
are v^^y = v sin^ cos 0 and Vn^z = AnO sin 0 , where 
angles 0 and (f characterize the incident direction of 
the wave packet [see Fig. 3(b)]. We denote the GH 
and IF shifts at the left (right) interface as A^^^ 

and A^^\ respectively. During multiple reflections 
in region II, the GH and IF shifts are accumulated, 
and induce average anomalous velocities: Vay = 
H- Aqjj) cos (f + {A^p + / ( 2 At) and 

^Gh) ^ ^ (^/F ^ ^?f) ^]/( 2 At). 

At = d/{v cosO) represents the propagating time be¬ 
tween two subsequent reflections. Therefore, the normal 
and anomalous velocity result in the total velocity: 

'^t,y{z) = '^n,y{z) +'^ 0 ,^( 2 ;) • 

We set the parameters E = 100 meV, V = 150 meV, 
h = 10/im, d = 50 nm, = vr = 1-2 x 10^ m/s, 
= V = 10^ m/s, = Vr = 0.8 x 10^ m/s, and 
velocities in other directions are all set as 10 ^ m /m 
Fig. i c) shows <p dependence of the total anomalous 
velocities induced by the spatial shifts. The anomalous 
velocities are valley-dependent, which implies that the 
total velocities also depend on valley index. Eventually, 
the different velocities of the two valleys lead to macro¬ 
scopic separation in real space, which is experimentally 
detectable. Fig. [^d) shows the opposite position shift 
at the bottom of region II for valley A and valley B. GH 
shift cannot induce the position shift, because GH shift 
always lie in the propagating plane. In contrast, IF shift 
can induce the position shifts (/im order) for two valleys 
(SN and QN)P. 

Identification of WSMs. —The experimental verifica¬ 
tion of WSMs has not been justified mainly due to 
the scant of efficient detection methocP^®^. The di¬ 
rect ARPES measurements of energy dispersions are cur¬ 
rently scant due to the constraint of magnetic properties 
of WSMs. Here, we suggest that the IF effect can be 
used as an experimental identification of WSMs. It has 
been shown that the topological IF shift splits the inci¬ 
dent wave packets into opposite directions respect to the 
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FIG. 4: Valley density distribution shift caused by IF effect, 
(a), (b) and (c) show the relative density distribution of valley 
A, valley B, and both, respectively, after considering GH and 
IF effect. The color bar on the right represents the relative 
density pjpmax of Weyl fermions. Here, p {pmax) represents 
the absolute (maximum) density of Weyl fermions at the bot¬ 
tom. The dashed lines indicate the location of the maximum 
intensity of density. Red (Green) box in Fig. (c) guides eyes 
to the region where only valley B (A) exist (pure valley po¬ 
larization). (d) Proposed setup for generating valley current. 
The devices Ti (red) and T 2 (green) are two terminals, which 
can extracted valley current out. 


valley index. Thus, the splitting of wave packets on the 
bottom of region II in Fig. [^a) can serve as a hallmark of 
WSMs. Furthermore, even considering an incident wave 
packet with finite angle range, this exotic splitting can 
also exist. 

Pure valley polarization. — Let us consider an TEM in¬ 
jector with incident angle range 6 {0 ^ [Oc — 5/2^ ^c+^/2], 
and 0 G [0c “ <^/2,0c + ^/2]) at the top of Region II, 
and study the valley density distribution at the bottom 
[see Fig. id)]. We calculate the valley density distri¬ 
bution with parameters E = 100 meV, V = 150 meV, 
h = 10 /im, d = SO nm, 6 = 3°^ Oc = 32.5°, and 0c = 90°. 
Without considering anomalous velocities, the density 
distributions for valley A and B are maximized at the 
center of y direction, and thus not distinguishable^. In 


contrast, after considering the anomalous velocities. Fig. 

(b) and (c) show the density distributions for valley 
A, valley B, and valley A & B, respectively. The IF effect 
induces opposite shifts in y-direction for valley A and val¬ 
ley B. Since the Weyl fermions from valleys A and B are 
well separated in space of micrometer order, pure valley 
current can be generated in the green (red) regions [see 
Fig. B c)]. This schematic set-up in Fig. [^d) can be 
utilized to generate pure valley current. 

Deteetion of Berry eurvature. —The Berry curvature— 
a gauge invariant quantity—should be detectable in 
experimenfP. However, up to date, there lacks an exper¬ 
imental feasible method to measure the Berry curvature 
in real materials. The topological IF effect provides a 
new way to measure the Berry curvature. For a system 
with inversion symmetry, the Berry curvature is an even 
function of wave vector k, i .e^ f2(k) = f2(-kP, which 
is usually satisfied in WSMs^^J®. Considering a wave 
packet propagating in x-y plane with kz = 0, the IF shift 
can be expressed as A/i? = 2 dk^ ^y(k), or 


^?/(^) ~ o 


1 9 A/i7(k) 


dkr 


( 7 ) 


To measure the Berry curvature one need to collect 
the IF shift in 2 : direction Ajf{E,0) as a function of 
energy E and incidence angle 0. The Aif{E^0) can be 
transformed into Aif{^xi and its derivative gives out 
the Berry curvature f^y(k). In the same way, the IF 
shift in other directions can be used to obtain ^^^^(k) and 
0.(k). 

Summary. —^We obtained analytical expressions of the 
GH shift and IF shift at the interface of WSMs. We 
demonstrate that the IF shift is valley-dependent, and 
can be attributed to the topological nature of the sys¬ 
tem. This IF shift can lead to valley-dependent anoma¬ 
lous velocity, which is experimentally detectable. Finally, 
we discuss three applications of the topological IF shift 
including characterization of WSMs, fabrication of high 
efficient valleytronic devices, and detection the Berry cur¬ 
vature. 
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Appendix A: Supplementary Information for “Topological Imbert-Fedorov shift in Weyl semimetals” 


Qing-Dong Jiang, Hua Jiang, Haiwen Liu, Qing-feng Sun, and X. C. Xie 

1. Detail derivation of the GH and IF shifts. 


To calculate the spatial shifts of the wave packet in ^ directions, one need to know the central positions of the 
incident and the reflected wave packets at the interface (x = 0). By expanding the phases a{ky,kz) and (j)r{ky^kz) to 

the first order around (ky, kz), the integrals of 'ipg^ and give that oc e ^ky/‘^ g-(^=F 2 g^j^d 


i^''g± OC e 


dky^l dky) ^ky!^ ^ ^ ^ ^ ^ ^ ^ Z ^ 


where the ± subscript corresponds to the first and the second 
component of the spinor, respectively. Then one can identify the center (maximum) of the incident and reflected wave 
packets in real space. For the incident wave packet, the two spinor components are centered at where 


' f£ = 

^2 

< 

; ; 

^in _ 
Z^ — 

^2 


(Al) 


ky - ky ,kz - kz 

For the reflected wave packet, the two spinor components are centered at {y±, where 


y± = 


^re 


dky ^r{ky, kz) T 2 dky^(^y^ 

~ at ^z) T 2 at 


ky=ky,kz=kz 

ky=ky ,kz=kz 

Considering above results, we can obtain the spatial shifts for two spinor components in y, z directions, which are 


(A2) 


= -^Mky, h) T h) 


= z'^ — z: 


-^re -^in 
± 


(A3) 


As is shown in Eq. (A3), in order to calculate the spatial shift at the interface (x=0), we need to know the reflection 
phase (j)r(ky^kz) and a{ky, kz). aiky^kz) represents the incident angle of the wave packet, and can be easily obtained 
once we know the incident wave vector. Next, we elaborate on deriving the reflection phase 0^. We consider the 
total reflection case, i.e., the reflection probability |rp = 1. Thus, the wave function must be evanescent in the region 
X > 0[Fig. 2(a)]. We can calculate the reflection coefficient by matching the wave function '0*^(k, r) + r) at 

X = 0 to the evanescent wave. The continuity of wave function gives the reflection amplitude r = where (pr = 

2 tan~^( _^ -a- 7r/2. In this expression, ^ = E-v+tlk, > where k = ^ {hv’ykyf + - {E - Vflhv'^. 


Substitute the expression of a and pr into Eq.(A3), and we can obtain the average shifts in y, z directions. 


II. IF shift for nonlinear energy dispersion. 


In the main text, we have shown that the Berry curvature is the origin of the IF shift for WSMs, which possess 
linear energy band dispersions. Here, we show that the Berry curvature is still the origin of the IF shift for another 
system without linear energy dispersions. To clarify this idea, we consider a system with a model Hamiltonian 


n' = 


Akl 

B{kx + ikz) 


B{k^ -ik^) 
-Akl 


(A4) 


where A and B are two parameters. The eigenenergy of the Hamiltonian is E = B'^{kl + fcf) + A^ky. Thus this 

Hamiltonian obviously does not have 3D linearized energy dispersions. One can apply a gate voltage V to this system 
at the region x ^ 0; therefore, an interface appears at x = 0. 

We now consider that a wave packet incidents on the interface from the region x < 0. To calculate the IF 
shift at the interface, we follow the same procedure in the main text. Analogously, we first construct the incident 
wave packet in Gaussian profile with center momentum (kx^ky^kz). Still we assume the wave packet as '0*’^(r) = 
f^oo dkydkz f{ky — ky) f{kz — kz) '0*^(k, r), where /{kg —kg) has the same definition as in the main text. Similarly, 
the reflected wave packet can be constructed as (r) = d kyd k^ f{ky - ky) f{k^ - k^) r). However, 
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since the Hamiltonian 1-L' here is different from that in the main text, one need to solve the wave functions r) 

and r) for . The Hamiltonian equation is l-L'ijj = where pj = (V^i, '^ 2 )^ is a two components spinor. Solve 

this Hamiltonian equation and use continuity condition of the wave function at x = 0, then we can get the solutions: 


V^-(k,r) 


1 

\/i + C" 


^-ia{ky,kz)l2 
^ ^Oiiky ^kz) j 


ikxX-\-ikyy-\-ikzZ 
^ 5 


and 


V>-(k,r) 


r 


_i ^iot{ky,kz)l2 
i ^-ia{ky,kz)l2 


^—ikxX+ikyU+ikzZ 


(A5) 


(A6) 


In 


/ E-ak?, 


the above expressions, ~ ~ {^^zY a = tan ^{kzl^x\ C = ^ ~ 

reffection amplitude, and (pr = 2 tan~^{ ) with and n = ^{E — Vy — {AkP^Y — [Bk^Y jB. 

Substitute (pr and a into Eq.(A3), and one can obtain the spatial shifts in y and 2 :-directions. If the incident wave 


packet is confined in x-y plane, then the out-of-plane spatial shift in ^-direction corresponds to the IF shift and reads: 

kr 


Aif = 


A Ekl 


(A7) 


Alternatively, one can also obtain the IF shift Ajf from semiclassical equation for wave packet (Eq.(7) in the main 
text). In order to apply the semiclassical equation, one need calculate the Berry curvature for the Hamiltonian 7^', 
where the superscript ± corresponds to the conductance band and valence band, respectively. The Berry curvature of 
this system is ky^. Substitute the Berry curvature expression into the Eq.(ll) 

in the main text, and one can obtain the analytic result of the IF shift Ajf = , showing the consistence with 


Eq.(A7). The fact that the IF shift Ajf can be obtained from Eq.(7) in the main text provides another convincing 
argument that IF shift arises from the Berry curvature of the system. 


HI. Calculations of anomalous velocities and position shifts for two valleys. 






FIG. 5: (a) The anomalous velocities induced by the GH shift are valley-independent. (b)The anomalous velocities induced by 
the IF shift are valley-independent. 


In the main text, we have demonstrated that the GH shift does not depend on valley index, whereas the IF shift 
depends on valley index. Although both of the GH shift and IF shift can bring anomalous velocities to the normal 
velocities, only the anomalous velocity induced by IF shift is valley-dependent. There are both GH and IF shifts at 
the left interface {x = — |) and the right interface (x = |). If the system has the mirror symmetry about the x = 0 
plane, the GH shifts are in the same direction at the left and right interfaces. But IF shifts are in opposite direction at 
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these two interfaces, leading to the cancellation of the IF shifts after two subsequent reflections. Therefore, in order to 
observe the valley-dependent IF shift, the mirror symmetry about x = 0 plane need to be broken. So we consider that 
the z direction velocities in the three regions are unequal, with v\ = vl^ and = vr. For simplicity, the x 

and y directions velocities are still with ~ Figure SI shows the anomalous velocities induced 

by the GH shift and IF shift. The anomalous velocities are calculated with the same parameters as in the section 
Anomalous velocities induced by IF effect in the main text. As one can see, the anomalous velocity induced by the 
GH shift is valley-independent; by contrast, the anomalous velocity induced by the IF shift is valley-dependent. The 
green curves and red curves correspond to valley A and valley B, respectively. The thick curves and dashed curves 
correspond to the anomalous velocities in y-direction and z-direction, respectively. 

Given the incident direction (6>, 0) of the incident wave packet, we can calculate the final position (x, y, z) at bottom 
of the region II. The wave packet needs time Tq = hjvn^z fo reach to the bottom of the region II without considering 
anomalous velocity. However, if we take the anomalous velocity into account, the wave packet needs time Ti = hjvt^z 
to reach to the bottom of the region H. Thus the anomalous velocities induced position shift is 


AT — ^t,y ^ Ff ^n,y ^ Fq. 


(A8) 


Put the total velocities of the valley A and valley B into Eq.(A8), we can get the final position shifts for the valley A 
and valley B, respectively [see Fig. 3(d) in the main text]. 


IV. Calculations of valley density distribution induced by valley-dependent anomalous velocities. 
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FIG. 6: Valley pattern before considering the IF shift. The figures (a) and (b) show the valley pattern of valley A and B, 
respectively, before considering the anomalous velocities induced by the spatial shifts. The color bar on the right represents 
the relative density of Weyl particles. Here, the relative density has the same meaning as that in the main text. 


Since the incident wave packets are likely to be reflected repeatedly at two interfaces, one only needs to consider the 
total reflected cases since the partially reflected wave packets would eventually diminish after being reflected multiple 
times. A collimated incident wave packet usually has certain angle range, thus, we assume 6> G [Oc — d/2^0c F d/2] 
and 0 G [(l)c — d/2^(j)cF d/2]^ where Oc and (j)c are the angles of the center of the collimated incident wave packet and 
the d is the angle spread range. In the numerical calculations, we take Oc = 32.5° and (pc = 90°. In this case, the 
anomalous velocity Va^y induced by the GH effect is negative, and the anomalous velocity induced by the IF shift is 
relatively large. Hence, the total anomalous velocity has a remarkable valley-dependent effect. 

For a wave packet of Weyl fermions propagating in the Region H, the total velocities in z directions are Vt^y and 
Vt^zi respectively. Due to the valley-dependent anomalous velocities, the total velocities are also valley-dependent, 
which can lead to macroscopic shifts of valleys. Because of the reflection at the two interfaces (Fig. 3(B)), the velocity 
in X direction in the Region H is a periodic function of At, i.e., {t) = {t + 2At). In the first period, the 

velocity {t) reads 


(() = 


t€ [0,f )U[5f‘,2Ai] 




t G 


[ At 

2 




(A9) 




















9 


The height of this structure is h, thus the wave packet needs time to = hjvt^z to reach to the bottom of the Region 11. 
Hence, one can get the final position coordinate (x^y^z) of the wave packet, where x = (itu^^(t), y = dtvt^y^ 

and z = —h. As one can see, a single incident direction (6>, 0) of the wave packet corresponds to a single point (x, y, z) 
at the bottom of Region 11. We further assume that the injector supplies an equal strength of incident beams in the 
angle range 0 G [32.5° — J/2,32.5° + 8/2] and (j) G [90° — ^/2,90° + 8/2]. Thus, there will be a pattern (density 
distribution) at the bottom of Region II [see Fig. 4(d) in the main text]. The Figure S2 (a) and (b) show the density 
distribution of valley A and valley B, respectively, before considering the IF shift. From the Figure SI (a) and (b) we 
can find that the density distribution of valley A and valley B are located at the same place. Therefore, the Weyl 
fermions from valley A and valley B are mixed together. However, if the anomalous velocity is taken into account, 
the distribution of valley A and valley B will shift to opposite direction leading to the valley separation [as shown in 
Fig. 4(a), (b) in the main text]. 
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